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We demonstrate the first solid-state spin-wave optical quantum memory with on-demand read¬ 
out. Using the full atomic frequency comb scheme in a Pr^+:Y 2 Si 05 crystal, we store weak coherent 
pulses at the single-photon level with a signal to noise ratio > 10. Narrow-band spectral filtering 
based on spectral hole burning in a second Pr®^:Y 2 Si 05 crystal is used to filter out the excess noise 
created by control pulses to reach an unconditional noise level of (2.0 ± 0.3) x 10“® photons per 
pulse. We also report spin-wave storage of photonic time-bin qubits with conditional fidelities higher 
than a measure and prepare strategy, demonstrating that the spin-wave memory operates in the 
quantum regime. This makes our device the hrst demonstration of a quantum memory for time-bin 
qubits, with on demand read-out of the stored quantum information. These results represent an 
important step for the use of solid-state quantum memories in scalable quantum networks. 

PACS numbers: 03.67.Hk,42.50.Gy,42.50.Md 


Photonic quantum memories are essential in quan¬ 
tum information science (QIS) where they are used as 
quantum interfaces between flying and stationary qubits. 
They enable the synchronization of probabilistic quan¬ 
tum processes e.g. in quantum communication [l|, 0 and 
computing Q. The implementation of quantnm memo¬ 
ries (QMs) for light requires strong interactions between 
individual photons and matter. This can be achieved by 
placing individual quantum systems (e.g. single atoms) 
in high finesse cavities Q or by using ensembles of atoms, 
where the photons are mapped onto collective atomic ex¬ 
citations. Atomic systems are natural candidates as QMs 
[Ml, but solid state systems offer interesting perspec¬ 
tives for scalability and integration into existing technol¬ 
ogy UMl- 

Rare-earth ion doped solids are promising candidates 
for high performance solid state QMs since they have 
excellent coherence properties at cryogenic temperatures 
m. They also exhibit large static inhomogeneous broad¬ 
ening of the optical transitions which can be tailored and 
used as a resource for various storage protocols, e.g. en¬ 
abling temporally 23| and spectrally 2^ multiplexed 
quantum memories. Recent experimental progress in¬ 
cludes qubit storage 15, H- 271, hig hly efficient quantum 
storage of weak coherent states [1^ , storage of entangled 

, entanglement between 


[l6j |. with a storage time limited by the coherence of ex¬ 
cited states. Another solution is to transfer the optical 


17, 18 


and single photons 

two crystals [i^ and quantum teleportation (30|. 

Yet, nonclassical states have so far only been stored 
as collective optical atomic excitations with fixed stor¬ 
age times [Illillii. While this may provide a use¬ 
ful resource if combined with massive multiplexing and 
deterministic quantum light sources 2^, the ability to 
read-out the stored state on-demand is essential for ap¬ 
plications where the quantum memory is used as a syn¬ 
chronizing device. On-demand read-out can be achieved 
by actively controlling the optical collective excitations 


spin-waves), using strong control pulses 1,3ll-3,311 . This 
gives access to much longer storage times [3J, • Op¬ 


erating a solid state spin-wave memory in the quantum 
regime has so far remained elusive, because of insufficient 
signal-to-noise ratio (SNR) at the single-photon level 

Here, using the full atomic frequency comb (AFC) pro¬ 
tocol in a Pr^+:Y 2 Si 05 Crystal we store and retrieve weak 
coherent pulses on-demand with SNR > 10 for single¬ 
photon-level input. Using a narrowband filter based on 
spectral hole burning in a second crystal, we achieve an 
unconditional noise floor of (2.0 ± 0.3) x 10“^ photons 
per pulse. The use of spectral holes as narrowband filters 
has been already demonstrated in storage schemes oper¬ 
ating in a classical fashion [sMi, but it is here exploited 
for the first time to enter the quantum regime. Finally, 
we demonstrate storage and retrieval conditional fideli¬ 
ties (i.e., assuming that a photon was reemitted) higher 
than classical memories for time-bin qubits at the single¬ 
photon level, taking into account the Poissonian statistics 
and the finite efficiency of the memory. These results rep¬ 
resent the first demonstration of a solid state spin-wave 
quantum memory, enabling on-demand read-out of the 
stored qubits. They also provide the first example of a 
spin-wave quantum memory for time-bin qubits (for any 
system), an essential resource in quantum communica¬ 
tion and processing [i^ . 

The AFC technique [H, is based on spectral tailor¬ 
ing of an inhoinogeously broadened absorption line into 
a comb-shaped structure with periodicity A. The in¬ 
put pulses resonant with the comb are mapped onto a 
collective optical atomic excitation. After an initial de¬ 
phasing, the excitations rephase at a time 1/A giving rise 
to a forward collective emission 15, 231. Before the co- 
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herent emission two strong control pulses are applied to 
transfer the excitation to and from a long-lived ground 
state to achieve the spin-wave storage of the input pulses. 
The full AFC scheme requires ions with at least three 
ground states, one being used as auxiliary state for op¬ 
tical pumping The spin-wave storage efficiency is 
given by rjsw = Vafc x Vt ^ where r]AFC is the ef¬ 
ficiency of the storage at the excited state and depends 
on the optical depth and comb finesse [1^, rjc accounts 
for the decoherence during the ground state storage, and 
rjT is the transfer efficiency of the control pulses. 

The realization of the full AFC scheme in the single¬ 
photon regime is challenging as the strong control pulses 
create noise which may dominate the weak signal re¬ 
trieved from the memory. Two main mechanisms are 
responsible for this noise, i) Spatial leakage from the con¬ 
trol mode into the input mode due to scattering from the 
optical surfaces and ii) interaction of the control pulses 
with residual population in the spin storage state due to 
imperfect optical pumping. The latter includes collective 
effects, such as free-induction decay (FID), or incoherent 
fluorescent emission. Note that four-wave mixing is not 
a dominant source of noise in our system, in contrast to 
Raman memories in atomic vapors 4^ (see ^ for de¬ 
tails). To reduce the noise we employ spatial, temporal 
and spectral filtering. The spectral hltering is challenging 
in Pr3+:Y2Si05, as the input and control frequencies are 
separated by only 10.2 MHz (see Fig. [Ijb)). As a narrow- 
band spectral filter we use a second Pr^+iY^SiO .tscry stal 


where we prepare spectral holes of variable width [3f 
The experimental arrangement and the relevant en¬ 
ergy level scheme of Pr^’*' at 606 nm are shown in Fig. 
[TJ The main laser beam at 606 nm (Toptica TA-SHG 
pro) is split into three to be used as input mode, fil¬ 
ter preparation, and lastly for control pulses and mem¬ 
ory preparation. They all pass through acousto-optical 
modulators (AOMs) in double-pass conhguration, driven 
by an arbitrary waveform generator (Signadyne), to cre¬ 
ate the necessary pulse sequences. The beams are then 
carried with polarization-maintaining single-mode opti¬ 
cal fibers to another optical table where the cryostat 
is located. The maximum available optical powers are 
about 20 mW, 3.5 mW and 150 for control, filter 
preparation, and input modes, respectively, measured in 
front of the cryostat. The frequency of the 606 nm laser 
is stabilized by Pound-Drever-Hall technique to a home¬ 
made temperature controlled Fabry-Perot cavity housed 
in a vacuum chamber. The input light is linearly polar¬ 
ized close to the optical D 2 axis to maximize the interac¬ 
tion with the Pr^+ ions. The measured optical depth 
of the Pr^+ transition at 606 nm is about 7 for both 
memory and filter crystal. In both cases, the inhomo¬ 
geneous linewidth is about 6 GHz. After the storage, the 
retrieved signal passes through different diffraction order 
modes (-1st and 1st) of two consecutive AOMs, acting 
as temporal gate before passing through the filter crys- 
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FIG. 1: (a) Quantum memory setup. The memory (MC) 
and filter (FC) crystals are located inside a liquid-free cooler 
(Oxford V14) operating at a temperature of 2.5 K. They 
are both 3 mm long and doped with a Pr®"*" concentration 
of 0.05%. The control and input beams are steered towards 
the memory with an angle of ~ 1.5°, leading to a exctinc- 
tion ratio of 10~®. The beam diameters on the crystal are 
280 pm and 90 pm for strong and input modes, respectively. 
The weak coherent states are prepared by attenuating bright 
pulses with variable neutral density filters (NDFs). A portion 
of the input beam is picked up before the NDF and sent to 
a photodiode (PDl) for the calibration of the mean photon 
number per pulse. A mechanical shutter is used to protect the 
SPD during the memory and filter preparation. HWP: half¬ 
wave plate; AOM: acousto-optical modulator; DG: diffraction 
grating; SPD: single-photon detector. The dashed red beam 
indicates the filter preparation mode, (b) Hyperfine splitting 
of the first sublevels of the ground ®H 4 and the excited ^D 2 
manifold of Pr^”*" in Y 2 Si 05 . 


tab A diffraction grating (DG) is then used to filter the 
noise not resonant with the crystal. The retrieved signal 
is coupled with 60% efficiency into a single-mode fiber 
for connection to the single photon detector (SPD, Pico- 
Quant r SPAD-20, detection efficiency rjd = 60%, dark 
count rate ^ 10 Hz). The total transmission of the input 
beam from the cryostat to the SPD is about 13 %. 

We tailor the AFG using optical pumping techniques 
as described in 33, li^. We isolate a single class of 
and create a 3.5 MHz wide AFG with A = 


44, 4£ 


atoms 

200 kHz on the l/2g — 3/2e transition within a 14 MHz 
wide transparency window (see ^ for details and for a 
comb example). During the preparation of the mem¬ 
ory, the population removed from the comb is stored 
in the auxiliary 5/2g state, while the 3/2g is emptied. 
To further remove unwanted residual population in the 
3/2g state, we apply an extra series of 100 pulses on 
the 3/2e — 3/2g transition after the comb preparation. 
We then start the single-photon-level storage measure¬ 
ments. Weak gaussian input pulses with full width at 
half maximum (FWHM) duration of 430 ns and mean 
photon number are mapped on the AFC and trans¬ 
ferred to spin-waves thanks to a strong control pulse. 
The control pulses have a Gaussian temporal profile with 
FWHM of 700 ns and are spectrally chirped by 5 MHz 
about the 3/2e — 3/2g transition. For each comb prepa¬ 
ration, 1000 storage trials are perfomed with a repetition 
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FIG. 2: (a). Time histograms of the retrieved photons mea¬ 

sured for different input photon numbers when a transparency 
window 2 MHz wide is prepared in the filter crystal. The 
input {fiin ~ 0.9) and the control pulses, as measured in 
photon counting and from a reference photodiode (PD2 in 
Fig. m, respectively, are also displayed. The chosen 0.7 /rs 
wide detection window is indicated by the dashed lines about 
the three-level echo; it includes ~ 80 % of the counts in the 
full echo mode. (b). Signal-to-noise ratio (SNR) as a func¬ 
tion of the number of input photons for different filter width. 
Circles: 2 MHz; Squares: 14 MHz. The error bars (smaller 
than the data points) are evaluated with Poissonian statis¬ 
tics. The black dotted line indicates the limit of detection 
SNR = 1. The dashed lines are linear fits of the experimental 
data. Panel (c). Decay of the SNR as a function of the spin- 
wave storage time Ts with average photon number /r™ = 1. 
From the fit with a Gaussian profile, the spin inhomogeneous 
broadening = (26 ± 1) kHz can be extrapolated. By com¬ 
paring the SNR that we measure at a storage time Ts = 7.8 /rs 
with the extrapolation at Ts = 0 /rs, we can evaluate the de¬ 
coherence term r/c to be about 75 %. 


rate of ^ 7 kHz. The full cycle has a period of 700 ms, 
including memory preparation and light storage. It is 
synchronized with the cryostat cycle to reduce the effect 
of vibrations. The sequence is then repeated 500 to 1000 
times to acumulate sufficient statistics. 

Fig. HJa) shows the time histograms of the retrieved 
photons with different /ri„. For this measurement, the 
crystal filter has a hole width of 2 MHz. The spin-wave 
storage time is T 5 = 7.8 fis, leading to a total stor¬ 
age time of Ts = 1/A + Ts = 12.8/rs. From the trace 
with = 0 , we estimate an unconditional noise floor 
of (2.0 ± 0.3) X 10“^ photons per pulse at the memory 
crystal. For fj,in = 1.15, we measure SNR = 16.3 ± 2.4. 
The linear scaling of the echo SNR with respect to in¬ 
creasing fiin is shown with blue circles in Fig. HKb). 
Typical values of efficiencies are tjafc = (5.6 ± 0.3) % 
and rjsw = ( 2.8 ± 0 . 1 ) %, from which we deduce r]T = 


(81.7 ± 2.6) % (assuming rjc 75% [S^). A convenient 
figure of merit taking into account the noise and effi¬ 
ciency is given by the minimum necessary to detect 
a spin-wave echo with SNR = 1, called fii. From the 
linear fit of the experimental data of Fig. [21(b), we find 
= 0.069 ± 0.002. We then vary by changing Ts 
with fiin = 1, as shown in Fig. He). The decay in the 
SNR is compatible with a spin inhomogeneous broaden¬ 
ing of 7m = (26 ± 1) kHzjSimilar to previous measure¬ 
ments with bright pulses 31, . We still observe SNR 


= 4.5 ±0.4 for Ts = 18.8 fis (r^ = 23.8/rs). 

We can estimate the contribution of the filter crystal in 
the suppression of the noise by preparing a wider trans¬ 
parency window. For a filter width of 14 MHz (squares 
in Fig. Hb)) there is no filtering at the control fre¬ 
quency, and we observe an increase of the noise floor 
to (2.3 ± 0.6) X 10“^ together with a slightly higher re¬ 
trieval efficiency {rjsw = (2.9 ± 0.2) %), which results in 
/ii values up to about 1. When the filter crystal is by¬ 
passed, the noise floor raises to (0.23 ± 0.01), indicating 
that the inhomogeneously broadened absorption profile 
of Pr^+ also contributes to partially filter the noise [ 2 ^ . 
Nonetheless, for this set of measurements, we were able to 
achieve higher storage efficiency, i.e. rjsw = (5.3±0.5) %, 
leading to a limited increase of the to about 4. 

For applications in QIS, it is crucial that the optical 
memory preserves the coherence of the stored qubits. We 
take advantage of the intrinsic temporal multimodality of 
the AFC protocol to demonstrate the phase preservation 
in the spin-wave storage of time-bin qubits. This type of 
encoding is widely used in quantum communication as it 
is robust against decoherence in optical fibers [s^. The 
time-bin qubits are expressed as \tpin) = ci|e)±C 2 e®^“|l), 
where |e) {\l)) represents a qubit in the early (late) time- 
bin, Aa is their relative phase, and cf ± C 2 = 1. In order 
to store time-bin qubits, the duration of the input pulses 
is reduced (from 430 to 260 ns), leading to a reduction of 
Vsw to about 2.2% and to an increase of the /ii up to 
flip = 0.11 ± 0.01. We start by evaluating the fidelity of 
the states |e) and |^), located at the poles of the Bloch 
sphere, Fg and F), by storing only the early and the late 
qubits. We obtain average fidelity values for the poles 
ranging from Fei = 85 % to 98 % for photon number per 
qubit, fiq, going from 0.6 to 5.9 [41 1. 

We then store superposition states located on the equa¬ 
tor of the Bloch sphere. We use the memory itself to an¬ 
alyze the retrieved qubits [s^ applying two partial write 
pulses as depicted in Fig. 121(a). This method provides a 
convenient way of analyzing time-bin qubits, but has the 
drawback of reducing the storage efficiency. As a matter 
of fact, in order to insert two write pulses, their duration 
needs to be reduced, which decreases their efficiency. Fig. 
|2Kb) reports examples of interference fringes for fiq = 1.5. 
From sinusoidal fits, we obtain a raw mean visibility of 

_= (72.5 ± 1.3) %. The fidelity of the process is calcu¬ 
lated from the visibility, asF=(l±F)/2 [^. Finally, 
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FIG. 3: (a) Pulse sequence to measure the time-bin qubit 
coherence. We apply two partial write pulses with a relative 
phase A/3 in order to split each pulse into two temporally sep¬ 
arated echoes. If the delay between the time-bins |e) and |/} 
is equal to the time difference between the two write pulses, 
we can overlap the late echo of the early bin, \el), with the 
early echo of the late bin, |Ze). The output of the memory 
(occurring after the single read pulse r) has three time-bins, 
{|ee), \el) -I- \le), |//)}. An interference will occur in the central 
time-bin if the coherence is preserved. The input pulses, lo¬ 
cated at |e) and |/) have a relative phase difference of Aq. (b) 
Interference fringes obtained integrating over the central out¬ 
put time-bin (in this case Atd ~ 0.5 fis) as a function of the 
relative phase difference A/3 for fiq = 1.5. Circles: Aa = 90 °, 
V = (71.6±6.8)%; Squares: Aq = 135°, % = (73.4±3.5)%. 


we obtain a total conditional fidelity per retrieved qubit 
Ft = ^Fei ± _, where Fei {F^ _) is the average fi¬ 

delity over the poles (equator) basis. The obtained values 
are reported in Fig. |4]for different ^q. We observe that 
the fidelity decreases with To explain this behavior, 
we fit our data with a simple model taking into account 
the decrease of SNR with fig and the reduced efficiency 
due to the double write protocol [^. The good agree¬ 
ment between the simple model and the data provides 
evidence that the decrease of fidelity is only due to the 
noise created by the control pulses, and not to a loss of 
coherence. 


In order to infer the quantum nature of our memory, 
the total fidelity is compared with the highest fidelity 
achievable with a measure-and-prepare approach (solid 
curve in Fig. n taking into account the Poissonian statis¬ 
tics of the input states and the finite memory efficiency 
(2.2%) 0, [HI, [ 2 ^. Using this criteria, the experimen¬ 
tal data are higher than the limit for a classical mem¬ 
ory by more than one standard deviation for most fiq 
investigated. With the raw data, the memory is in the 
quantum regime for fiq > 0.96. When correcting for the 
loss of efhciency in the analysis (which could be achieved 
by analyzing the qubits with an external interferometer), 
the model predicts that the quantum regime would be 


reached for ^q > 0.25 (see 0). 

The very low noise probability and the ability to ob¬ 
tain <C 1 opens prospects for the spin-wave storage 
of single-photons, as required for many applications in 
QIS. In that case, the probability to have a photon be¬ 
fore the memory (i.e including all optical loss between 
source and memory) needs to be higher than to en¬ 
ter the quantum regime. In the current experiment, the 
storage efficiency is limited by the available comb optical 
depth in our 3 mm long crystal, by the limited trans¬ 
fer efhciency {rj^ = (67 ± 4) %) and by technical issues 
(cryostat vibrations, laser linewidth, see 0). Note that 
much higher efficiencies (for storage in the excited state) 
have been obtained in Pr^+:Y 2 Si 05 Using longer crystal 
(69 %) 0 or impedance matched cavities (58 %) (0 . 
Longer storage times can also be achieved with dynam¬ 
ical decouplin g te chniques to counteract decoherence in 
the spin state 34, ^ £7, ^|. 



FIG. 4: Total fidelity vs input photon number per qubit, 

The light orange squares are the data points with an error 
bar of 1 standard deviation. The orange dotted line is a fit to 
the data points using Eq. (5) in f4il| . with the corresponding 
shaded area being the 1 standard deviation of the error in this 
fit. The solid blue (dashed green) line is the classical limit 
obtained by a measure and prepare strategy for our memory 
efficiency of rjsw = 2.2 % (r/sw = 100 %) when testing the 
memory with weak coherent states [^. The dash-dotted line 
is the classical limit for testing the memory with a single¬ 
photon Fock state (F = 2/3). 


In conclusion, we demonstrated the spin-wave storage 
and on-demand retrieval of weak coherent states at the 
single-photon level in a solid state memory based on a 
Pr^+:Y 2 Si 05 Crystal. This is the first demonstration that 
solid state spin-wave optical memories can operate in the 
quantum regime, overcoming a strong limitation for AFC 
QMs. We achieved a SNR higher than 10 for single¬ 
photon level input pulses, high enough to enable the stor¬ 
age of single photons. Finally, we confirmed the quantum 
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nature of our memory by storing time-bin qubits encoded 
in weak coherent states and demonstrating conditional fi¬ 
delities for the retrieved qubits higher than what is pos¬ 
sible with classical memories. Our device thus represents 
the first spin-wave memory for photonic time-bin qubits. 
These results open the door for long-lived storage and 
on-demand readout of non-classical light states in solid 
state devices and represent an important step in view of 
using solid state quantum memories in scalable quantum 
architectures. 
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Supplemental Material 

In the supplemental material we report details on 
the atomic frequency comb preparation (section L), the 
sources of noise (section II.), the chosen strategy to sup¬ 
press the technical noise (section III.), the characteriza¬ 
tion of the filter cristal (section IV.), and, finally, we de¬ 
scribe a theoretical model for the fidelity of the spin-wave 
storage of time-bin qubits (section V.). 


Atomic Frequency Comb Preparation 


To prepare the atomic frequency comb (AFC) we fol¬ 
low the approach already described in 33|, l43|. The 
frequency of the strong beam is firstly scanned by 14 
MHz in 100 ms to create a wide transparency window 
(also referred to as pit) in the Pr^+ absorption line by 
spectral hole burning (curve a in Fig. 0. Afterwards 
a 2 MHz-wide sweep is operated outside the pit to pump 
the atoms back to the 1 /2g state. This creates a 3.5 MHz¬ 
wide absorbing feature corresponding to the l/2g — 3/2e 
transition inside the pit, but also populates the 3/2g 
state which in principle must be empty for later trans¬ 
fer of the input field (see Fig. 1(b) in the main text). 
Thus, a clean pulse with duration of 50 ms is applied 
on 3/2g — 3/2e transition to empty this ground state. 
The clean pulse also contributes to suppress the tran¬ 
sitions of different classes of Pr^+ ions which might be 
addressed by the preparation pulses [i^. Then, a se¬ 
quence of single-frequency hole-burning pulses is applied 
on the single class absorption feature resonant with the 
l/2g — 3/2e transition, each time changing the frequency 
by a fixed amount, A. This sequence burns periodi¬ 
cally spaced holes in the absorbing feature corresponding 
to the l/2g — 3/2e transition and anti-holes at the fre¬ 
quency of the 3/2g — 3/2e transition, so a short burst of 
clean pulses is applied to maintain the 3/2g state empty. 
For the AFC preparation it is crucial to have a third 
ground state (5/2g) to use as auxiliary state where to 
store the excess of population resulting from the optical 
pumping. An example of the resulting comb structure 
for A = 200 kHz is shown in Fig. [Sj curve b. With 
the described procedure we are able to tailor absorp¬ 
tion peaks as narrow as (43 ± 3) kHz, with a comb op¬ 
tical depth d = 4.5 ± 0.1 and background optical depth 
do = 0.75 ± 0.04. It is worth noting that the limited dy¬ 
namical range of the detector might fix a higher bound 
for the detected peak amplitude and thus lead to an over¬ 
estimated full width at half maximum. 

When the memory is characterized with classical 
pulses, we find typical efficiency values of tjafc = (10.3± 
1.4)% and rjsw = (4.4 ± 0.1)%. The prediction for 
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FIG. 5: Example of transparency window (curve a) and comb 
trace prepared with A = 200 kHz (curve b). The inset shows 
a magnihcation of the comb corresponding to the l/2g — 3/2e 
transition. 

Gaussian AFC peaks is 'qAFC = , where 

d = d/F is the effective comb optical depth [23|. The 
measured value of tjafc is compatible with the predic¬ 
tion, rjAFCAh = (12.2 ±1.6) %, for the measured values of 
d = 4.5±0T, F = 4.7±0.4 and do = 0.75±0.04. For the 
same optical depth of our comb but assuming do = 0, the 
model predicts a maximal efficiency of r]AFC,th = 26.5 % 
for a finesse of F = 4. Note that an even higher effi¬ 
ciency, riAFC,th = 34.8%, is expected for square shaped 
comb peak and a finesse of F = 3.2 (). In our case, the 
achievable contrast is limited by technical reasons, such 
as vibrations of the cryo-cooler and finite linewidth of 
the laser. When measured in photon counting mode, the 
efficiency is further reduced by the coincidence window 
taken (containing 80% of the total counts), by loss in 
the narrow-band hlter ('^ 10%) and by instabilities due 
to long integration times (^10% loss). 

Sources of noise 

The full AFC protocol for light storage implies the use 
of strong transfer pulses which are sources of different 
kind of noise. First of all there might be spatial leakage 
from the preparation/control mode into the input mode, 
essentially due to reflections and scatting from the opti¬ 
cal surfaces. Moreover, the control pulses might interact 
with the residual population left in the 3/2g state, due 
to imperfect cleaning (see section ,). This interaction can 
give rise to coherent (e.g. free-induction decay) or inco¬ 
herent emission. Another potential source of noise in this 
system may arise from four-wave mixing, where the first 


control pulse creates spin-waves via an off-resonant Ra¬ 
man interaction and the second control pulse reads them 
out as noisy photons in the same mode as the memory 
output. This has been identified as the main source of 
noise in Raman type quantum memories using atomic 
vapors [i^. However, this type of noise does not seem 
to be, at this stage, an important source of noise in our 
spin-wave memory. This is evidenced by the following 
noise measurements. First, we have measured that the 
two control pulses contribute equally to the noise in the 
echo window. This can be explained by the fact that they 
interact with the same number of ions giving two equally 
sized incoherent fluoresecence fields (given the time scale 
of the storage time compared to the excited state lifetime 
Ti = 164 [m). If four-wave mixing was not negligible, 

the noise would be suppressed by more than a factor 2 by 
switching off one of the two control pulses. Furthermore, 
the noise arising from four-wave mixing would decay for 
increasing storage times with the same behavior of the 
spin-wave echo, i.e. following the inhomogeneous spin 
broadening (see Fig. 2(c) of the main text). As a mat¬ 
ter of fact, the noise floor in our case is constant within 
error over the whole range of storage times investigated 
(see Fig. [6|). 
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FIG. 6: Noise floor in the echo window as a function of the 
time between the two control pulses, F. The error bars are 
calculated taking into account Poissonian statistics. 

Several reasons could explain this difference. First, the 
non-collective emission from the ^D 2 excited state mainly 
happens through intermediate crystal-field manifold and 
involves the coupling with the phonons of the YSO ma¬ 
trix. As a matter of fact, the branching ratio of the 
spontaneous emission through the ^D 2 —?>^H 4 transition 
is very low (Ri 3% ) jHlj. Furthermore, among the tran¬ 
sitions within the ^H 4 —>■^02 that could be off-resonantly 
excited by the first control pulse, many of them are 
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characterized by relative oscillator strengths one order 
of magnitude lower than that employed to transfer the 
excitation to the ground state. During the AFC prepara¬ 
tion, most of the ions are collected in the auxiliary 5/2g 
state but, due to the relative oscillator strength of the 
transitions starting from it, they can only be involved 
in the almost closed cycle transition 5/2g o 5/2e Q. 
Those characterized by a non-negligible probability for 
off-resonant excitation either are suppressed by the state 
cleaning that we operate in the memory preparation (it 
is the case of the transitions starting from the 3/2g state) 
or are tailored in the comb structure (transitions having 
the l/2g as initial state). In the latter case, the second 
control pulse would give rise to a coherent emission in a 
well separated temporal mode with respect of the spin- 
wave echo. 


Filtering 

To achieve a low noise floor, (2.0 ± 0.3) x 10“^, we 
use spatial, temporal and spectral filtering. As for the 
spatial filtering, the control and input beams are sent to 
the memory with a small angle, leading to an extinction 
of 10“®. Temporal filtering is achieved thanks to a tem¬ 
poral gate implemented with two AOMs. It allows us 
to block the leakage of the strong control beams in the 
signal mode. This is important not to blind the single¬ 
photon detectors, but also to avoid burning a spectral 
hole in the crystal filter. A diffraction grating is em¬ 
ployed to decrease the noise originated from incoherent 
fluorescence. The narrow-band spectral filter is prepared 
simultaneously with the AFC. The frequency of the filter 
mode (shown in dashed red line in Fig. 1(a)) is scanned 
by 1.2 MHz around the input light frequency for 100 ms 
which results in a transparency window of around 2 MHz 
due to the power broadening effect. In order to tempo¬ 
rally discriminate the stored and retrieved weak pulses 
from the portion of the FID leaking in the echo mode 
and happening at short time scales (a few fis), the op¬ 
timization of the excited state storage efficiency at long 
storage times is crucial. By carefully optimizing the stor¬ 
age times, in both the excited and ground states, and the 
pulse durations a compromise can be reached between the 
noise level and a reasonable retrieval efficiency. 


Filter crystal characterization 

We provide in this section details about the noise sup¬ 
pression realized by the filter crystal. The main advan¬ 
tage of using a second crystal is that it provides a dynam¬ 
ical filter, whose center frequency and bandwidth can be 
easily tuned, as it is obtained by spectral hole burning. 

Fig. [7|a) represents the noise floor (noise counts per 
pulse at the memory crystal) as evaluated from noise 
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FIG. 7: Noise floor in the echo window (a) and fii (b) value 
as a function of the spectral hole width in the Hlter crystal. 
Filled symbols: hole centered at the input pulse frequency; 
open symbols: transparency window coincident with that pre¬ 
pared in the memory crystal (see curve b in Fig. (Sjl. (c) Time 
histogram of the retrieved photons for different input pho¬ 
ton numbers, fj,in, and a filter width of 14 MHz. The input 
(fim = 0.9) and the control pulses (measured with a reference 
photodetector) are also displayed. 


measurements where we burn in the filter crystal spectral 
holes of different width centered about the input pulse 
frequency (solid symbols). The noise in the echo win¬ 
dow increases by one order of magnitude when the filter 
width increases from 0.8 to 12 MHz. The open symbol 
refers to the extreme case of a 14 MHz-wide transparency 
window coincident with the one prepared in the memory 
crystal to host the AFC (see curve a in Fig. [5]). Here the 
noise floor is further doubled because the transparency 
window allows the control pulse frequency to pass by and 
reach the detector. The increase in the noise affects con¬ 
sequently the lii value as shown in Fig. EKb). Neverthe¬ 
less, this value also depends on the echo counts, thus the 
narrowest filter width does not necessarily correspond to 
the best case. As a matter of fact, the best /ii value is 
measured for a filter width which allow the major por¬ 
tion of the echo signal to pass by but suppresses efficiently 
the noise at the control pulse frequency, i.e. 2 MHz. The 
value of jjLi also depends on the echo efficiency and is 
thus more affected by power fluctuations than the noise 
floor. We believe this is the reason why its increase in 
the case of totally open transparency window is more 
pronounced than that of the noise floor (compare filled 
and open symbols in panel (a) and (b) of Fig. [7]). Panel 
(c) of Fig. [7] shows the time histogram of the retrieved 
photons for different input photon numbers in the case of 
a 14 MHz-wide pit prepared in the filter crystal. Finally 












we evaluate the extinction ratio of onr 2 MHz-wide spec¬ 
tral filter by switching off the AFC and sending input 
pulses through the transparency window in the memory 
crystal firstly resonant with the filter center, then at the 
control pulse frequency. By comparing the total counts 
reaching the detector in the two cases, we estimate the 
extinction ratio to be about 750. 


Modeling the Fidelity 


Fig. 4 from the main document shows the measured 
total fidelity as a function of the mean photon number 
per qubit, fiq, with a corresponding fit using the model 
we develop in this section. The aim of the model is to 
show that the reduction of the fidelity as the photons 
per qubit are reduced is entirely due to the noise of the 
memory and not due to a loss of coherence caused by 
e.g. phase noise in the laser. The goal of our model is 
to express the total fidelity as a function of the same fii 
parameter. 

Firstly, let us look at the output of the memory for the 
two differing measurements, namely the measurement of 
the |e) or |/) qubit (i.e. the poles of the Bloch sphere) 
and the measurement of the coherence between |e) and 
|/) (i.e the equators of the Bloch sphere). The former 


(latter) are schematically depicted in Fig. 8(a) (8(b) |. 


For a given integration width, we sum over both E 
and L for the poles measurement, while for the coherence 
measurement we integrate over the window labelled E- 
L-l-L-E (see figure [8(b)| ) . 


Consider the example depicted in Fig. 8(a) where the 
|e) qubit has been stored and then recalled. We label the 
integrated counts in the window about E as S -I- N and 
the integrated counts about the window L as N, where S 
is the signal (not including the noise) and N is the noise. 
Then, the fidelity is the ratio of the counts in E to the 
total counts over E and L, i.e. 


T. = 


S-tN 


S/N-hl SNR-hi 
S-hN -h N “ S/N-h2 “ SNR + 2’ 


( 1 ) 


where SNR is the signal to noise ratio. Provided the 
efficiency and noise are the same for storing the \l) qubit, 
Ti= Te = J-ei- Finally, we can express the above fidelity 
as a function of the /ri of this measurement (see main 
text for definition) by noting that SNR = p,g/^ip, where 
/rip is the /ri characteristic of the poles measurement. 
The expression for the fidelity becomes 


^el = 


/Tg ~\~ /Tip 

/ig -\- ‘2fi\p 


( 2 ) 


and we measure /rip = 0.11 ± 0.01. 

Now let us consider Fig. |8(b)[ which depicts an ar¬ 
bitrary qubit |e) -I- \l) being recalled, and the relative 
phases are such that the maximum interference is ob¬ 
served. Firstly, the integrated counts corresponding to 



(a) 



(b) 


FIG. 8: (a) Schematic example of the memory output for 
storage of an |e) qubit, (b) Schematic example of the memory 
output for storage of an |e) -I- \ l) where the relative phases are 
such that the maximum interference is observed. 


the window at the time labelled E-E are S/4a-|-N where 
the signal S is the same as defined earlier. Here, the signal 
S is reduced by a factor of 4, this is because the number 
of photons per qubit are distributed equally within the 
input time-bins (factor of 2) and the echo is distributed 
equally between the output time-bins (factor of 2). The 
factor a reduction in the signal is because we use differ¬ 
ent bandwidth pulses for the double write process which 
gives a reduced efficiency in the echo. The amount of sig¬ 
nal in the window E-L -|- L-E is 4 times that of the signal 
in the L-L window provided the interference is fully co¬ 
herent (note the same is true for the E-E window).The 
resulting counts are then S/a -|- N. It is worth noting 
that if we instead used an unbalanced Mach Zehnder in¬ 
terferometer, we would have no reduction in the signal 
and a = 1. 


The fidelity of the coherence measurement is ob¬ 
tained from the visibility of the interference fringe by 
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T = (1 + F)/2 [ 3 ^ where the visibility V is defined as 


max — min S/a + N — N 

max + min S /a + N + N 

S SNR 

S + 2Na “ SNR + 2a 


(3) 


where max (min) refers to the maximum (minimum) in¬ 
tegrated counts of the fringe, the max being S/a-l-N and 
min being N . Note that we have V in terms of the SNR 
of the poles measurement, so the resulting fidelity can be 
written in terms of /rip, i.e. 


= i _l_ 1 

2 2 /Tg “t“ 2a//ip 


(4) 


where the subscript -l-(—) refers to the qubit |e) ± |Z) 
(|e)±i|^)), and we assume = J'- = We measure 
a by comparing the fii in the L-L window of figure [s/b)] 
to /rip and we get a = 2.5 ± 0.6. 


concluding our simple model for the fidelity. Taking the 
measured flip, we fit the data in Fig. [^with equation 
[5] having a as a free parameter. We get a = 2.5 ± 0.3 
which agrees excellently with the measured value stated 
earlier. We therefore can conclude that the reduction 
of the fidelity is due to the noise on the output of the 
memory and not phase noise. For fig > 0.96 the fidelity 
surpasses the classical bound, remaining in the quantum 
regime. Also in Fig. [9] we plot the case for a = 1 for 
flip = 0.11 and flip = 0.07 (the case shown in Fig. 2 of 
the main text). The fidelity exceeds the classical bound 
for photon numbers much less that 1 for these cases. 

For completeness, we include table U which shows all 
the measured fidelities for every fig tested. Also, we in¬ 
clude examples of the output of the memory for early and 
late inputs (Fig. [TUT all and for the input |e) -I- \l) where 
the phase of the double-write process is adjusted to show 
maximum and minimum interference in the E-L + L-E 
window. 



FIG. 9: Total fidelity vs input photon number per qubit. The 
light orange squares are the data points with an error bar of 
1 standard deviation. The light orange dotted line is a fit to 
the data points using eq. [5] with flip hxed at 0.11 and a — 
2.5 ± 0.3. The corresponding shaded area is the 1 standard 
deviation of the error in this fit. The solid blue (dashed green) 
line is the classical limit obtained by a measure and prepare 
strategy and a memory efficiency of t/sle = 2.2 % (r/sLE = 
100 %) when testing the memory with weak-coherent states 
am. We see that for fig > 0.96 the fidelity is above the 
classical benchmark. We show also the case of a = 1 and 
flip ~ 0.11, i.e. the case where our detection of the coherence 
is replaced by an unbalanced Mach Zehnder interferometer for 
example. Furthermore, we show the case of a = 1 and flip = 
0.07, i.e. the fii measured for the optimised pulse bandwidth 
case. The dash-dotted line is the classical limit for testing the 
memory with a single-photon Fock state {J- = 2/3). 


Now, the total fidelity is given by J^t = \^ei + f-Tq_, 

which leads to 


Tt 


1 

3 


flq T f^lp 

fig ‘2fllp 



flq ‘ley.flip _ 


(5) 


flq 

J=.i 

J-+ 

Tt 

Tc 

5.9 

(97.9 ± 1.5)% 

(97.2 ± 1.7)% 

(97.4 ± 1.2)% 

(93.0 ± 0.1)% 

3.2 

(96.9 ± 2.2)% 

(93.3 ± 1.2)% 

(94.5 ±1.1)% 

(90.1 ±0.1)% 

1.5 

(93.5 ± 2.9)% 

(86.2 ± 1.6)% 

(88.6 ± 1.4)% 

(86.2 ±0.1)% 

1.1 

(93.3 ±3.1)% 

(85.8 ± 1.5)% 

(88.3 ± 1.4)% 

(84.4 ± 0.1)% 

0.6 

(84.9 ± 4.2)% 

(72.9 ± 1.3)% 

(76.9 ± 1.6)% 

(81.0 ±0.1)% 


TABLE I: Fidelities T^i, Tj, _and Tt for each /t, tested. Tc 

is the classical benchmark for testing the device with weak 
coherent states which depends on the efficiency of the memory 
and the input photon number used i,m- The stated errors 
are one standard deviation. 


* Present address: Department of Physics, Clarendon Lab¬ 
oratory, University of Oxford, Oxford 0X1 3PU, United 
Kingdom 

^ Electronic address: margherita.mazzera@icfo.es 

[1] L.-M. Duan, M. D. Lukin, J. 1. Cirac, and P. Zoller, Na¬ 
ture 414 , 413 (2001). 

[2] N. Sangouard, C. Simon, H. de Riedmatten, and N. 
Gisin, Rev. Mod. Phys. 83 , 33 (2011). 

[3] E. Knill, R. Laflamme, and G. J. Milburn, Nature 409, 
46 (2001). 

[4] H. P. Specht et al, Nature 473 , 190 (2011). 

[5] C. W. Chou et al, Nature 438, 828 (2005). 

[6] T. Chaneliere et al., Nature 438 , 833 (2005). 

[7] K. S. Choi, H. Deng, J. Laurat, and H. J. Kimble, Nature 
452 , 67 (2008). 

[8] H. Zhang et al, Nat Photon 5, 628 (2011). 

[9] M. Lettner et al, Phys. Rev. Lett. 106 , 210503 (2011). 

[10] S. Riedl et al, Phys. Rev. A 85, 022318 (2012). 

[11] A. Nicolas et al, Nat Photon 8, 234 (2014). 

[12] B. Julsgaard et al.. Nature 432, 482 (2004). 

[13] M. Hosseini et al., Nat Phys 7 , 794 (2011). 























10 



FIG. 10: (a) An example of the detected counts at the output 
of the memory for the early (light orange histogram) and late 
(dark blue histogram) input qubits with ^.q = 3. The vertical 
dashed lines correspond to the 0.5 /rs integration windows for 
the calculation of the fidelity, (b) An example of detected 
counts at the output of the memory for an input qubit |e) + |Z) 
with /Tg = 5.9. The dark blue (light orange) histogram shows 
the case when the relative phase of the double-write process 
is such that the interference is maximized (minimized) in the 
integration window (indicated by the vertical dashed line). 
The pulse feature on the right corresponds to the L-L window 
of Fig. |8(b)| while the E-E window is not detected here as the 
temporal gate is closed to negate the effects of the control 
pulses. 


[14] M. Sprague et ai, Nat Photon 8, 287 (2014). 

[15] H. de Riedmatten et ai, Nature 456, 773 (2008). 

[16] M. P. Hedges, J. J. Longdell, Y. Li, and M. J. Sellars, 
Nature 465, 1052 (2010). 

[17] C. Clausen et al., Nature 469, 508 (2011). 

[18] E. Saglamyurek et al, Nature 469, 512 (2011). 

[19] E. Togan et al., Nature 466, 730 (2010). 

[20] H. Bernien et al., Nature 497, 86 (2013). 

[21] D. G. England et ai, Phys. Rev. Lett. 114, 053602 
(2015). 

[22] R. M. Macfarlane, J. Lumin. 100, 1 (2002). 

[23] M. Afzelius, C. Simon, H. de Riedmatten, and N. Gisin, 
Phys. Rev. A 79, 052329 (2009). 

[24] N. Sinclair et al., Phys. Rev. Lett. 113, 053603 (2014). 


[25] M. Giindogan et al., Phys. Rev. Lett. 108, 190504 (2012). 

[26] C. Clausen, F. Bussieres, M. Afzelius, and N. Gisin, Phys. 
Rev. Lett. 108 , 190503 (2012). 

[27] Z.-Q. Zhou et ai, Phys. Rev. Lett. 108, 190505 (2012). 

[28] D. Rielander et ai, Phys. Rev. Lett. 112, 040504 (2014). 

[29] 1. Usmani et al, Nat Photon 6, 234 (2012). 

[30] F. Bussieres et al, Nat Photon 8, 775 (2014). 

[31] M. Afzelius et al., Phys. Rev. Lett. 104 , 040503 (2010). 

[32] M. Giindogan et al. New Journal of Physics 15, 045012 
(2013). 

[33] N. Timoney et al, Phys. Rev. A 88, 022324 (2013). 

[34] J. J. Longdell, E. Fraval, M. J. Sellars, and N. B. Manson, 
Phys. Rev. Lett. 95 , 063601 (2005). 

[35] G. Heinze, C. Hubrich, and T. Halfmann, Phys. Rev. 
Lett. Ill, 033601 (2013). 

[36] H. Zhang et al.. Applied Physics Letters 100, 131102 

( 2012 ). 

[37] S. E. Beavan, M. P. Hedges, and M. J. Sellars, Phys. Rev. 
Lett. 109 , 093603 (2012). 

[38] S. E. Beavan, E. A. Goldschmidt, and M. J. Sellars, J. 
Opt. Soc. Am. B 30 , 1173 (2013). 

[39] 1. Marcikic et al.. Nature 421, 509 (2003). 

[40] P. Humphreys et al., Phys. Rev. Lett. Ill, 150501 
(2013). 

[41] See Supplemental Material, which includes Refs. [4y-[^. 

[42] P. S. Michelberger et al.. New Journal of Physics 17 , 
043006 (2015). 

[43] N. Maring et al.. New Journal of Physics 16, 113021 
(2014). 

[44] M. Nilsson et al, Phys. Rev. B 70, 214116 (2004). 

[45] O. Guillot-Noel et al., Phys. Rev. B 79, 155119 (2009). 

[46] M. Sabooni, Q. Li, S. Kroll, and L. Rippe, Phys. Rev. 
Lett. 110 , 133604 (2013). 

[47] M. F. Pascual-Winter, R.-C. Tongning, T. Chaneliere, 
and J.-L. Le Gouet, Phys. Rev. B 86, 184301 (2012). 

[48] M. Lovric, D. Suter, A. Ferrier, and P. Goldner, Phys. 
Rev. Lett. Ill, 020503 (2013). 

[49] O. Guillot-Noel et al., Phys. Rev. B 75 , 205110 (2007). 

[50] M. Bonarota, J. Ruggiero, J.-L. Le Gouet, and T. 
Chaneliere, Phys. Rev. A 81, 033803 (2010). 

[51] D. L. McAuslan, J. J. Longdell, and M. J. Sellars, Phys. 
Rev. A 80 , 062307 (2009). 

















